Non-local effects in effective medium response of nano-layered meta-materials 
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We analyze electromagnetic modes in mult i- layered nano-composites and demonstrate that the 
response of a majority of realistic layered structures is strongly affected by the non-local effects 
originating from strong field oscillations across the system, and is not described by conventional 
effective-medium theories. We develop the analytical description of the relevant phenomena and 
confirm our results with numerical solutions of Maxwell equations. Finally, we use the developed 
formalism to demonstrate that multi-layered plasmonic nanostructures support high-index volume 
modes, confined to deep subwavelength areas, opening a wide class of applications in nanoscale light 
management. 



Nanolayered composites have been recently proposed 
to serve as negative index systems, super- and hyper- 
lenses, photonic funnels, and other nanophotonic 
structureo 1 ! 2 ! 3 ! 4 ! 5 ' 6 ! 7 ' 8 ! 9 ' 10 ' 11 , The typical thickness of an 
individual layer in these "artificial" (meta-) materials is 
of the order of lOnm. Since this size is much smaller 
than optical (or IR) wavelength, it is commonly as- 
sumed that the properties of the multilayered compos- 
ites are well-described by the effective medium theory 
(EMT) 13 i 14 . In this Letter, we analyze the modes of re- 
alistic multi-layered structures and show that the con- 
ventional EMT fails to adequately describe these modes 
due to the metamaterial analog of spatial dispersion - 
strong variation of the field on the scale of a single layer. 
We derive a non-local correction to EMT, bridging the 
gap between metamaterial- and photonic crystal-regimes 
of multi-layered media, and use numerical solutions of 
Maxwell equations to verify our results. Finally, we use 
the developed technique to identify volume metamaterial 
modes confined to nanoscale areas. 

While the formalism developed below is applicable to 
the composites with arbitrary values of permittivities op- 
erating at different frequency ranges (UV, visible, IR, 
THz), here we illustrate our approach on the optical re- 
sponse of a two-component plasmonic nanolayered com- 
posite, which has been suggested for a variety of future 
beam-steering and imaging systems 3 -! 4 -! 1 - 1 -. The schematic 
geometry of such a structure, containing alternating lay- 
ers of materials with permittivities, t\, 62 and (average) 
thicknesses a\ and ci2 respectively, is shown in FigfTJ In 
the analytical results presented below, we mostly focus 
on the propagation of TM waves, which are responsible 
for plasmon-assisted phenomena; and only briefly discuss 
the implications for TE modes. It is straightforward to 
generalize the presented technique for mixed waves as 
well as for multi-component structures. In the selected 
geometry, x coordinate axis is perpendicular to the plane 
defined by layer interfaces, while y and z axes are paral- 
lel to this plane; the direction of z axis is chosen so that 



FIG. 1: Schematic geometry of a planar nanolayer-based 
meta-material, surrounded by two cladding layers 



that the electromagnetic waves propagate in x, z plane. 

The majority of realistic designs of layered nanoplas- 
monic structure s 3 ' 4 ' 11 ! 12 rely on the metamaterial regime, 
when the typical layer thickness is much smaller than the 
free-space wavelength A so that surface plasmon polari- 
tons propagating on different metal-dielectric interfaces 
are strongly coupled to each other. The optical prop- 
erties of the metamaterial structure in this strong cou- 
pling regime could be related to some effective permittiv- 
ities. The existence of these effective parameters is im- 
portant from both fundamental and applied standpoints. 
Thus, effective-medium description provides one with an 
insight into the physics behind the optical response of the 
structure; it also significantly reduces the computational 
efforts needed to simulate the electromagnetic response 
of the system, simultaneously increasing the accuracy of 
these simulations; finally, it can be easily utilized to pro- 
vide an important link between the properties of (easily 
fabricatable) systems with a few layers, with the ones of 
(more practical) macroscopic multi-layered structures. 

Apart from the wavelength, two more independent 
length-scales can be identified in the system - the one 
of the typical layer thickness a ~ 0,1,0,2, and the one of 
the typical field variation L. Since the introduction of 
effective permittivity e requires some kind of field av- 
eraging, independence of L and A yields to a fundamen- 
tal difference between the metamaterial 1 - 3 - and "conven- 
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tional" effective-medium 14 responses of nano-composites. 
As we show below, in nanoplasmonic layered structures 
L < A so that e cff will have non-local corrections. 

As any optical system, the multi-layered composite 
can be described by the behavior of its resonant (eigen) 
modes. Each such mode is characterized by its effective 
modal index, given by n e s = k z c/u> with k z and c being 
the modal wavevector and speed of light in the vacuum 
respectively. An arbitrary wave propagating through the 
system can be represented as a linear combination of dif- 
ferent modes. Note that when n 2 s > exceeds that of 
both cladding layers, the electromagnetic field of a mode 
is confined inside the layered structure, which behaves 
like a waveguide. 

To analyze the electromagnetism in the metamaterial, 
we numerically solve 3D Maxwell equations in the lay- 
ered geometry using the Transfer Matrix Method. In this 
technique, described in details in 1 ^, the field in each layer 
is represented as a combination of two (plane) waves hav- 
ing the same dependence in z direction and propagating 
in the opposite x directions, followed by the construction 
of a transfer matrix describing the collective response of 
the multi-layered structure. The modes of the metamate- 
rial are then related to the eigen values and eigenvectors 
of the transfer matrix. 

To understand the evolution of multilayered system 
between meta-material and effective-medium regimes, 
we used the transfer-matrix techniques to identify the 
modes of a 200-nm-thick layered composite with per- 
fectly conducting cladding layers - essentially represent- 
ing a waveguide with deep subwavelength crossection. 
This technique allows us to control the field variation in 
the direction perpendicular to the waveguide, and simul- 
taneously enforce the "metamaterial condition" 0.1,2 
A = 1.55 fim for all nanolayered structures in our work. 

We generated ~ 100 ensembles of nanocomposites with 
ei = —100 and e 2 ~ 2 (Au/SiC>2 composite). In each 
ensemble, we fixed the total thickness of the compos- 
ite h = 200nm, total concentration of metal, and ran- 
domly varied thickness of individual metal layers. The 
variation in layer thickness was about 10% of average 
thickness. The idea behind the ensemble generation is 
two-fold. First, we aim to understand the response of re- 
alistic multi-layered systems, where the total number of 
layers is relatively small; second, this approach gives us 
the opportunity to assess the tolerance of the composite 
properties with respect to fabrication defects. 

The profiles of several eigen modes are shown in Figj^l 
Note that the field across an individual layer is expo- 
nential rather than oscillatory in nature, so that high- 
precision arithmetic is required to find the accurate nu- 
merical solution of Maxwell equations. As the number 
of layers is increased (and correspondingly the thickness 
of an individual layer is decreased), the field distribu- 
tion in the system converges to the one of the mode in 
a waveguide with homogeneous core. Therefore, in this 
regime the behavior of a multilayered composite is essen- 
tially identical to the behavior of a uniaxial anisotropic 




FIG. 2: (Color online) TMM (solid lines) and EMT (dashed 
lines) calculations of E z field of TMi(a,b), and TM 2 (c,d) 
modes for metal (white) -dielectric (yellow) composites with 
Ni = 10 (a,c) and 20 (b,d) layers 



system with effective permittivity tensor e , given by 
< D a >= < Ep >, with Greek indices correspond- 
ing to Cartesian components and <> being the average 
over the multi-layer subwavelength area^. Due to axial 
symmetry, e cff is diagonal, its optical axis coincides with 
x, and ef = ef = ef. 

The dispersion relations of the TM and TE waves prop- 
agating in such a metamaterial 11 : 

<f _ kl kl + kl 

r 2 f eff ,eff ' 

yz x 

_ k 2 , + ky + k 2 
72" — lil ' \ > 

respectively (as noted above, k y — 0). 

We calculated waveguide modes for each composite in 
an ensemble. The results of our numerical solutions of 
Maxwell equations and their comparison to conventional 
EMT— with 

£ off = £ (0) _ (ai +Q 2 )eie 2 

02£l + a\€2 

fiff _ c (0) _ + Q2C2 , , 

e yz — e yz ~ . \ z ) 

are summarized in Fig|3] It is clearly seen that simi- 
lar to what has been shown for fiber geometry ir*£, the 
planar multi-layered composite supports highly-confined 
volume modes. It is also seen that while the response 
of all structures in a single ensemble is very alike, and 
therefore the introduction of effective permittivity is jus- 
tified, conventional EMT fails to describe the behavior of 
majority of practical nanolayered composites. A reason- 
able agreement is present only when the number of layers 
Ni = 2h/(a 1 + a 2 ) is very large. Note that the EMT does 
not work despite the fact that the condition «C A is 
met. 

The origin of this effect lies in a strong variation of the 
fields on the scale of a single layer, clearly visible in Fig. 2. 
Similar to the strong field variation on the subatomic 
scale that yields non-local corrections to permittivities 
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of homogeneous materials^, the scale separation L < A 
introduces non-locality into e off . Note that in contrast to 
the case of non-local response in homogeneous structures, 
the microscopic (layer-specific) field in meta-material can 
still be described by "local" ei j2 ; "effective" non-locality 
is present only in the effective permittivity. 

To find the non-local correction to the EMT, we start 
from the layered metal-dielectric structure where all 
metallic and all dielectric layers have the same thickness 
(a i may be still different from a 2 ). In this limit, the 
system essentially becomes a ID photonic crystal (PC). 
The dispersion of the modes of this case can be related 
to the eigen- values problem for two-layer transfer matrix, 
yielding^ 

s(k x [ai + a 2 ]) = cos(fc 1 a 1 ) cos(fc 2 a 2 ) 
— 7sm(/e 1 ai) sm(k 2 a 2 ) 



cos( 



(3) 



where the polarization-specific parameter 7 is given by 
1 / e 2 k\ , t\k 2 



7tm - 2 



e 1 k 2 e 2 fci 



7T E = 5(| + |) (4. 



and k\ 2 = ei, 2 ^ - k\. 

The '"conventional" EMT regime [Eqs.@] can be ob- 
tained from Eq.Q through the Taylor expansion up to 
the second order in |feiai| <C 1; \k 2 a 2 \ <C 1; \k x (ai + 
a 2 )\ <C 1 (see e.g.—). Expanding the PC dispersion equa- 
tion up to the next non-vanishing Taylor term yields se- 
ries of modes with dispersion given by Eq.([T]) and effec- 
tive permittivities 

M 



1 - 5 x (k,uj) 



,eff 



£ (0) 

tyz 



yz i-6 v ,(k, u ) 

where the nonlocal corrections are given by: 
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k 2 x {ei+e 2 f 
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(6) 



Note that since components of the wavevector are related 
to each other via Eq.{T]), the choice of k x and w/c as 
opposed to k z or k y [in Eq.©] is somewhat arbitrary 
and primarily depends on the geometry. Here we use 
k x = 2irj /h and cu as independent variables for j th mode. 

The agreement between the non-local EMT with re- 
sults of numerical solutions of Maxwell equations is 
shown in Fig[3J It is clearly seen that Eq.© perfectly 
describes the behavior of lower-order modes. The agree- 
ment tends to worsen for ?i c ff 3> I where \k ■ a\ > 1. To 
confirm that cladding regions and material absorption 
have weak effect on non-localities we have also simulated 
the modes of realistic Au-SiC>2 structures with vacuum 
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FIG. 3: (Color online) Effective refractive index of waveg- 
uide modes as a function of number of layers in the core 
region AT; = 2h/(a\ +02), calculated using transfer ma- 
trix method (dots), "conventional" EMT (dashed lines), and 
the "non-local EMT" derived in this work (solid lines), for 
h — 200nm-waveguide with perfectly conducting (a,b) and 
air (c,d) claddings; ai = 02/3 (a,c); a\ = 0,2 (b,d); top panels 
correspond to ei = — fOO; size of black bars correspond to 
standard deviation in n e ff as determined from our numerical 
simulations; bottom panels represent Au — Si02 composite 
with ei = e Au = -114.5 + ll.Olz 



cladding regions. The results of these simulations and 
their comparison to our non-local EMT are shown in 

FigEM). 

An important note is that "real" parameter behind 
the validity of effective medium response is \k ■ a\ <C 1. 
fn majority of all-dielectric nanostructures | ei.2 1 — 1 
or \k x c/uj\ 1, and this parameter is identical to the 
commonly- used criterion a A. For the high- index TM 
modes in metal-dielectric systems the spatial dispersion 
provides a significant correction to the quasi-static EMT 
results. The effective non-locality will be present for 
all-dielectric materials provided that |ex,a | 3> 1- Sim- 
ilar effect have been recently discovered for microwave 
nanowire structures^. 

To conclude, we have demonstrated that conventional 
EMT fails to adequately describe the optical proper- 
ties of multi-layered metal-dielectric metamaterials. We 
identified strong variation of the field to be the cause 
of this disagreement and derived an analytical correc- 
tion to incorporate non-local effects into EMT. We have 
also demonstrated that multi-layered structures support 
high-index modes confined to spatial areas as small as 
A/8. Our results, illustrated here for TM waves in two- 
component optical structures, can be used to design 
nano-guiding systems, and to provide an efficient link 
between the properties of realizable few-layer structures 
and their multi-layered macroscopic counterparts. The 
presented techniques are directly applicable to UV, IR, 
or THz metamaterials and can be generalized to mixed 
(HE, EH) waves and to multi-component systems involv- 
ing anisotropic materials using techniques of 15 . 
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